Partitioning of the linear photon momentum in multiphoton ionization 
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The balance of the hnear photon momentum in multiphoton ionization is studied experimentally. 
In the experiment argon and neon atoms are singly ionized by circularly polarized laser pulses with a 
wavelength of 800 nm and 1400 nm in the intensity range of 10^* — 10^^ W/cm'^. The photoelectrons 
are measured using velocity map imaging. We find that the photoelectrons carry linear momentum 
corresponding to the photons absorbed above the field free ionization threshold. Our finding has 
implications for concurrent models of the generation of terahertz radiation in filaments. 
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The photon momentum is usually eclipsed by other 
photon properties, such as energy or angular momen- 
tum. At high photon energies the photon momentum is 
manifest, e.g., in the weak non-dipole contributions to 
the photo-electron angular distributions [1]. Although 
the momentum of a single visible photon is exceedingly 
small, a large number of photons can give rise to macro- 
scopic effects of radiation pressure [2, 3]. In multiphoton 
ionization at relativistic light intensities radiation pres- 
sure is known to cause a suppression of recollision [4] and 
therefore the suppression of both non-sequential multiple 
ionization [5, 6], and high harmonic generation [7]. Re- 
cently, radiation pressure has been invoked as an expla- 
nation for the terahertz radiation emitted by laser gen- 
erated filaments [8-11]. However, the detailed balance 
of the quantized momentum transfer from the photon 
field to the electron-ion system in multi-photon ioniza- 
tion (MPI) is both theoretically and experimentally un- 
explored. 

We have studied the photo-electron momentum dis- 
tribution produced in circularly polarized, infrared, fem- 
tosecond laser pulses along the direction of the laser pulse 
propagation. We observe and quantify the effect of radia- 
tion pressure on isolated electrons following multi-photon 
ionization. We find that the asymptotic electron momen- 
tum in the direction of photon propagation corresponds 
to a fraction of the total momentum absorbed in multi- 
photon ionization. Specifically, neither the photons driv- 
ing the quiver motion of the free electron during the laser 
pulse nor the photons required to lift the electron from 
the bound state into the continuum contribute to the 
forward momentum of the free electron after the pulse. 
We will show that the former is clear from semi-classical 
physics. The physics underlying the latter has never been 
considered. 

Previous studies of the photon momentum transfer 
were performed at relativistic intensities and the long 
pulse regime. In these conditions the ionization poten- 



tial Ip is negligible and the photo-electron exits the focus 
in the presence of the high intensity pulse [12, 13]. In our 
experiment the photo-electron kinetic energy is compa- 
rable to Ip and the photoelectron remains in the focus for 
the duration of the pulse. These conditions are similar 
to most THz generation experiments [14]. 

In multi-photon ionization, the total energy absorbed 
by the photo-electron is given by the number of pho- 
tons [15] Nhu! = Ip + K + Up, where Ip is the ioniza- 
tion potential of the atom, K is the drift kinetic energy 
of the electron measured after the event and Up is the 
ponderomotive energy of the electron. In a circularly po- 
larized pulse Up = QeEo/ (2ma;^), where Eq and uj are 
the light electric field and frequency, respectively. Photo- 
electron distributions of MPI clearly reveal the quantized 
nature of the transfer of energy and angular momentum. 
Likewise, the transferred linear momentum must be con- 
served Nhk = pI +pI where z is the direction of light 
propagation and the superscripts denote the momenta of 
photo-electron, photo-ion, and laser field, respectively. 

We used circularly polarized laser pulses at 800 nm, 
15 fs (intensity FWHM) and 1400 nm, 70 fs to ionize 
Ar or Ne atoms. The energy of the 800 nm pulse was 
controlled using a half-wave plate and two germanium 
plates positioned at Brewster's angle. The energy of the 
1400 nm pulse was controlled using a neutral density 
filter. Our experiment spans the intensity range from 
10^** — 10""^^ W cm~^. At these intensities and wave- 
lengths, the photo-electron drift velocity is less than 2% 
of the velocity of light. 

We measure the photo-electron momentum using a 
velocity map imaging (VMI) spectrometer [16]. The 
spectrometer projects the 3-dimensional electron momen- 
tum distribution onto a 2-dimensional imaging detector 
(Fig. 1(b)). The detector plane contains the direction 
of laser propagation (z) and the axis of the gas jet (x). 
The laser pulses were focused into the super sonic gas jet 
with a local density of « 10^ cm~^. The velocity dis- 
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FIG. 1. (a) Photoionization of a neutral atom gives a net mo- 
mentum to the pliotoelectron pi = K/c and the ion pi = Ip/c 
in the direction of light propagation, (b) Measured electron 
momentum spectrum (in atomic units) of Ne in circularly po- 
larized light at 800 nm, 8 x fO" W/cm^. The laser pulse 
propagates along the z-axis. (c) Background image under the 
same conditions as in (b) reveals a distinct peak of low en- 
ergy electrons from laser produced Rydberg atoms that are 
ionized by the spectrometer DC field (indicated by the arrow), 
(d) The measured Ne photoelectron distribution (crosses) is 
compared to the Rydberg reference distribution (dots) and a 
reference Gaussian distribution centred at pz = 0. The centre 
of the Ne distribution has a net pz > 0. 



tribution of the atoms in our pulsed gas jet corresponds 
to a spread in the electron momentum distribution of 
< 1.5 X 10~^a.u. (atomic units - a.u.). 

The nominal spectrometer momentum resolution, de- 
termined by its pixel size, is 7 x lO^'^ a.u. for Ar and 
1 X 10~^ a.u. for Ne. However, multi-photon ion- 
ization creates a Gaussian wave packet in momentum 
space along the laser propagation k-direction [17] that 
is much broader than the nominal resolution. A typical 
wavepacket that we measure is shown as the crosses in 
Fig. 1(d). Taking the centroid of this distribution, we 
can determine a shift between momentum distributions 



obtained under differing conditions with an accuracy of 
~ 2 X 10"'^ a.u. Since the linear momentum of a single 
800 nm photon is 4 x lO^"' a.u., we are sensitive to a mo- 
mentum shift caused by the absorption of ~ 5 photons. 

To determine the absolute shift we use electrons ex- 
tracted from Rydberg states that are populated by multi- 
photon excitation [18, 19] as a reference. The lifetimes 
of these Rydberg states can easily exceed the dura- 
tion of our laser pulses. When the laser pulse passes, 
atoms left in highly excited Rydberg states can be ion- 
ized by the spectrometer's DC electric field, resulting 
in very narrowly distributed, low momentum electrons 
(see Fig. 1(c)). The Rydberg states are most efficiently 
populated in linear light. However, the water vapor in 
the background gas provides Rydberg electrons even in 
circular polarized light. We estimate that the electro- 
static fields of the spectrometer 0.6 kV/cm < Egtat < 
1.3 kV/cm can field ionize Rydberg states with n > 20. 

Rydberg electrons provide an unambiguous reference 
for the zero momentum. The Rydberg atom or molecule 
has absorbed approximately TV = Ip/{hw) photons. Since 
the multi-photon absorption and ionization steps are de- 
coupled in Rydberg atoms, the linear momentum Ip/c is 
transferred to the center of mass of the atom [19]. Hence, 
the origin of the Rydberg electron distribution coincides 
with zero velocity in the lab frame. By recording many 
Rydberg electron spectra for a range of laser intensities 
(shown in Fig. 2b) we determine the origin of the mo- 
mentum distribution with an accuracy of 1 x 10~'^ a.u. 

Our choice of circularly polarized light has two im- 
portant advantages over linear polarization. It avoids 
Coulomb focusing [20] in the electron momentum distri- 
butions and provides a simple measure for the laser in- 
tensity. In circular polarization the electron momentum 
distribution resembles a torus [21] whose radius is pro- 
portional to the peak of the laser field experienced by the 
neutral atoms. Since we are only interested in single ion- 
ization this is the appropriate intensity measure [22, 23]. 
We determine the light intensity with a precision of 8% 
[24], thereby directly measuring the net number of pho- 
tons absorbed for a typical electron at that intensity. 

For the data analysis we integrated the measured 2-D 
spectra over x to arrive at a 1-D distribution along the 
axis of laser propagation (the crosses in Fig. Id). We 
use a nonlinear fitting routine to fit the center, width, 
amplitude, and noise-floor of a Gaussian distribution to 
the measurements. 

In Fig. 2 the intensity dependence of the center of the 
distribution, dN/dpz, is shown for multi-photon ioniza- 
tion in 800 nm and 1400 nm of neon (circles) and argon 
(triangles, crosses). The intensity dependent momentum 
transfer has a different slope for the two wavelengths as is 
expected from the wavelength dependence of the asymp- 
totic electron kinetic energy K = q^E^/ (2rna;^) [25]. 
Generally, with increasing intensity the center of dN/dpz 
is increasingly shifted in the direction of propagation. 
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velope results in a net negative force as the pulse passes 
over the electron. 
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FIG. 2. (a) The center of the electron momentum distribution 
along z as a function of the average intensity at ionization for 
Ar and Ne at 800 nm and 1400 nm. The intensity at ionization 
is determined from the electron px distribution averaged over 
the ionization volume (space) and laser pulse length (time; 
see text for details). The theory lines are given by Eq. (2). 
(b) The zero position of the momentum distribution along the 
z axis as a function of laser intensity. The position is deter- 
mined by the center of the Rydberg electron spot indicated 
by the arrow in Fig. Ic. 



The upper limit of intensity arises because of saturation 
of single ionization, when a further increase of pulse en- 
ergy does not result in an increase in the intensity seen 
by the neutral atoms. The actual saturation intensity de- 
pends on the ionization potential and the pulse duration 
(70 fs at 1400 nm, and 15 fs at 800 nm). 

In the following we will analyze the transfer of 
photon momentum to a free electron in the clas- 
sical description of the light field [25], using the 
fact that the photon momentum can be described 
by the Lorcntz force [7]. A free electron is driven 
by the laser electric field (propagating along z) 
E(z, t) — Eog (z, t) [cos{ujt — kz) Cx + sn\{ujt — kz) ey\ 

with the pulse envelope g (z, t) = exp ^—^5^ (z — ctf'^ , 
T being the pulse length FWHM in intensity. With the 
magnetic field 3 = 6^ x E/c the equation of motion for 
the free electron is F = qE + (/v x B. At non-relativistic 
velocities the magnetic field and the slowly varying pulse 
envelope can be treated as perturbations on the electron's 
motion. To first order the electron's velocity v only de- 
pends on the electric field. The second order correction 
for the momentum component Pz includes the contribu- 
tion of the magnetic field, B, and the ponderomotive gra- 
dient VUp. The action of the magnetic field generates a 
net forward momentum. The gradient of the pulse en- 
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Note, that B is for a constant pulse envelope, whereas 
the second term —d^Up includes the pulse envelope. For 
intensities in the tunneling regime the atom is most likely 
to ionize around the peak of the envelope (t=0). Since 
the pulse is short we can set z = during the time the 
electron interacts with the laser. Using this simplifica- 
tion, Eq. 1 can be integrated to yield the net longitudinal 
momentum after the pulse has passed: 
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The first term represents the forward momentum due to 
the radiation pressure on the free electron. The second 
term is the effect of the gradient in Eq. 1. It removes the 
forward momentum associated with the ponderomotive 
energy. This is the reason Up does not appear in Eq. 2. 
When the pulse has vanished there is still a forward mo- 
mentum due to the electron's drift kinetic energy, K. 
We confirmed the validity of Eq. 2 using a numerical 
simulation to the complete Lorentz equation of motion 
including the pulse envelope. 

Eq. 2 agrees with the experimental data in Fig. 2. For 
comparison, in Fig. 2, we also plot Eq. 2 plus a constant 
offset Jp/c. The dashed line, which includes this offset, 
overestimates the measured p| well beyond our experi- 
mental uncertainty. Thus, we conclude that the portion 
of the longitudinal momentum /p/c, corresponding to the 
photons necessary to overcome the ionization energy, is 
not given to the photo-electron. It must be transferred to 
the centre of mass of the electron-ion system {p\ = Ip/c). 

The agreement between Eq. 2 and the measurement 
allows us to draw another conclusion. It shows that the 
ponderomotive energy does not transfer any net momen- 
tum to the electron. This is because the ponderomotive 
term in Eq. 1 removes that part of the longitudinal mo- 
mentum as the pulse passes the electron. The plasma- 
induced phase modulation of the pulse [26] returns the 
momentum — Up/c to the photon field. If the pon- 
deromotive energy were necessary, the slope in Fig. 2 
would be different by a factor of two - well outside the 
error bars. This conclusion means that linearly and cir- 
cularly polarized pulses will transfer very different longi- 
tudinal momentum to the photoelectrons. In linearly po- 
larized light most of the energy absorbed by the electron 
upon ionization is ponderomotive energy and is therefore 
returned to the field. In circularly polarized light, half 
the absorbed energy is in Up and half in K . The momen- 
tum associated with K remains with the photoelectron 

In Fig. 3 we compare the photoelectron longitudinal 
momentum calculated for circularly and linearly polar- 
ized light at two different wavelengths. In both cases 
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FIG. 3. Net photoelectron longitudinal momentum (pi) vs 
laser intensity calculated for linear and circularly polarized 
light. The curves are calculated from Eq. 2 using the kinetic 
energy averaged over one optical cycle. 



we have used Eq. 2 to calculate the longitudinal momen- 
tum. Circular and linear polarization have very different 
kinetic energy distributions and this affects the net longi- 
tudinal momentum of the photoelectrons when averaged 
over one optical cycle. In circularly polarized light, the 
ionization probability is independent of the phase of the 
field. In linear polarization, the ionization probability 
and the kinetic energy are strongly depend on the phase. 
An electron ionized at the phase ujti has drift kinetic en- 
ergy K — q^E^ sin^ ojti/ (2r7iw^) . The average kinetic en- 
ergy in linearly polarized light is K = J^^^^ K (ti) ^ dti 
where ti is the time of ionization and dN/dti = W [E (ij)] 
is the Yudin-Ivanov sub-cycle ionization rate [27]. In gen- 
eral K in linearly polarized light is much less than in 
circular polarization at a constant intensity. 

There are two important implications of our results. 
First, there has been a debate about the proper role of 
tunneling in multi-photon ionization [28, 29]. It is clear 
that tunneling - a DC phenomenon - provides no mech- 
anism to transfer linear momentum to the ion. Although 
tunneling is a tremendous computational tool, our exper- 
iment shows that it remains an approximation. To our 
knowledge, how the momentum Ip/c is shared between 
the electron and ion has never been addressed before in 
experiments, simulations, or models. 

Second, the photon momentum transfer to the elec- 
tron has been invoked to account for THz generation in 
filaments [8, 9, 11]. In the low fields typical of atmo- 
spheric gas breakdown, the momentum transfer will be 
very small for linearly polarized light. Our results sug- 
gest this can be augmented by using circular polarization 
and longer wavelength. 
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